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The time evolution of a charged point particle is governed by a second-order integro-differential 
equation that exhibits advanced effects, in which the particle responds to an external force before 
the force is applied. In this paper we give a simple physical argument that clarifies the origin and 
physical meaning of these advanced effects, and we compare ordinary electrodynamics with a toy 
qq , model of electrodynamics in which advanced effects do not occur. 
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Charged point particles in classical electrodynamics obey an equation of motion, known as the Lorentz-Dirac 

equation, that is third order in proper time, and therefore presents various difficulties [l|, 0, S 0] ■ in particular, the 

Lorenz-Dirac equation admits unphysical runaway solutions, in which the particle accelerates without bound even 

in the absence of an external driving force. The runaway solutions can be eliminated by replacing the Lorenz-Dirac 

i - i, equation with a second-order integro-differential equation, but this equation has problems of its own: it exhibits 

i-pH . advanced effects, in which the particle responds to an external force before the force has been applied. 

In this paper we give a simple physical argument that clarifies the origin and physical meaning of these advanced 

c/3 • effects, and we compare ordinary electrodynamics with a toy model of electrodynamics in which advanced effects do 

not occur. The toy model is closely analogous to ordinary electrodynamics; for example, there are fields that are 

i— H , analogous to the electric and magnetic fields, and these fields both mediate an interaction between charged particles 

and support freely propagating radiation. The toy model, however, does not exhibit the conceptual problems that 

Y^ plague ordinary electrodynamics; for example, the self-energy of a point particle does not diverge, and there are no 

. — , ■ advanced effects. Thus, the toy model shows that it is possible to construct a mathematically consistent theory of 

coupled particles and fields that is not subject to these pathologies. By performing a detailed comparison between 

^~. a the toy model and ordinary electrodynamics, we can gain insight into the conceptual foundations of electrodynamics, 

and into the meaning of the advanced effects. 

The paper is organized as follows. In section [Til we briefly review the Lorentz-Dirac equation of motion, the 

corresponding integro-differential equation of motion, and the causality problems associated with it. In section HTT1 

we consider a spatially extended charged particle, and give a simple physical argument that explains the origin of the 

advanced effects. In section ITVl we introduce the toy model of electrodynamics. In section [V] we consider a spatially 

extended particle in the toy model, and derive analogs to the various equations of motion for electrodynamics. 

(f) i The following notation is used in this paper. The function e(x) is the sign function, defined such that e(x) = 1 

if x > 0, e(x) = if x = 0, and e(x) = — 1 if x < 0. The function 9(x) is the step function, defined such that 

r — , 9(x) = 1 if x > 0, 6(x) = 1/2 if x = 0, and 9(x) — if x < 0. The metric tensor r/^ is defined such that 7700 = 1, 

7711 = 7/22 = 7733 = —1- The three- vector component of a four- vector is denoted by boldface type; for example, 

OO! x» = (x°,x). 
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Consider a point particle coupled to the electromagnetic field. We will let m and e denote the mass and charge of the 
particle, and we will let z^(t) denote the position of the particle at proper time r. Also, we will define v^ — dz^/dr 
to be the velocity of the particle and a^ = dv^ /dr to be its acceleration. The equation of motion for the particle is 

ma» = K» + K&t, (1) 

where K*t is the force exerted on the particle by the electromagnetic field and K^ xt is an arbitrary externally imposed 
force. For simplicity, we will assume that K^ xt depends only on the proper time r, and not on the position or velocity 
of the particle. The force K^ is given by 

K^r)=eF^(z(r))v v (T), (2) 



II. THE LORENTZ-DIRAC EQUATION 



where F 1 *" = d^A 1 ' — d v ' A^ is the electromagnetic field strength tensor and A^ is the vector potential. In the Lorentz 
gauge (d^A^ — 0), the vector potential satisfies the field equation 

DA^ = 4ttJ^, (3) 

where J^(x), the current density, is given by 

J»(x) = e I ' v»{T)5 iA) (x-z{T))dT. (4) 



Equations {THI} give a complete description of the coupled particle-field system. 

We can decompose each solution to the field equation ([3|) into the sum of an inhomogeneous solution, which describes 
the potential generated by the particle, and a homogeneous solution, which describes freely propagating radiation. It 
is useful to perform this decomposition in two different ways: 

A = A r + A in = A a + A out , (5) 

where A% and A£ are the retarded and advanced potentials generated by the particle, and A^ n and A^ ut describe 
incoming and outgoing radiation. The retarded and advanced potentials are given by 

A?(x)= [ D r (x-x').P l (x')d 4 x', A»(x)= ( D a (x-x')J»(x')d A x', (6) 



where D r (x) and D a (x) are the retarded and advanced Green functions for the inhomogeneous wave equation: 

D r (x) = 29(x°) S(x ■ x) = Ixl" 1 5(x° - |x|), D a {x) = 28(-x°) S(x ■ x) = Ix^ 1 5(x° + |x|). (7) 

Using the decompositions of the vector potential given in equation ([5]), we can express the field-strength tensor and 
the electromagnetic force as 

F»» = F,r + Ft: = K v + FZ , K» s = K? + K? n = K^ + K» ut , (8) 

where, for example, 

Fi? = d»A v r - VA?, K?(t) = eF^{z(r)) v v (t). (9) 

Physically, K n describes the force exerted on the particle by incoming radiation, and K r describes the self-force 
exerted on the particle by its own retarded field. In what follows, we will assume that there is no incoming radiation, 
so Kf n = and K^ f = K£. 

We can obtain an explicit expression for the self- force by combining equations (|4|), ([6]), and ([9]). To perform the 
calculation it is useful to express the self-force in the form K^ = Kit + Kt, where K± = (1/2) (K£ ± Kg). As was 
first shown by Dirac [i[ (see also Q, pages 187-9), the component Kt is well-defined and is given by 

Kt = mr (d M + (a • a)^), (10) 

where d M = da^/dr and tq = (2/3)(e 2 /m). The component Kt, however, is infinite. The reason for this can be traced 
to the fact that the retarded fields of a point particle diverge as we approach the particle; for a spatially extended 
particle, one can show that Kt is finite and is given by 

K% = -msa" + ■■■ (11) 

where ms is the self-energy of the particle and the dots indicate additional terms that vanish in the point particle 
limit. We derive an expression for the self-energy of a spherically symmetric particle in Appendix \^ but for now 
we simply note that the self-energy scales like 1/cr, where a is the particle size, and therefore diverges in the point 
particle limit. Thus, in order to obtain a finite expression for the self-force, let us assume that the particle is spatially 
extended, but small enough that Kt is well-approximated by equation (|10p and Kit is well-approximated by the first 
term of equation (fTTj) : 

K? = Kt + I<1 = mr (d M + (a • a)i/*) - msa*. (12) 

Substituting this result into equation (fl]) , we find that the equation of motion for the particle is 

ma^ = mT (d M + (a ■ a)v M ) - m,sa M + Kg xt . (13) 



Here m is the bare mass of the particle; that is, the mass that the particle would have if its charge were set to zero. 
Let us define a renormalized mass m R = m + ms and a renormalized time constant t r = (m/m R )TQ — (2/3)(e 2 /m R ). 
We can then express the equation of motion as 

a" - r R a» = TR (a ■ a)v" + {l/m R )K» xt . (14) 

This is the Lorentz-Dirac equation. We derived this equation by considering a small spatially extended particle, but it 
is well-defined in the point particle limit, provided we hold m R constant and allow m to diverge in order to compensate 
for the divergence of ms- 

It is instructive to consider the low-velocity limit of the Lorentz-Dirac equation: 

a - t r sl = (l/m R )K ext , (15) 

where a = da./dt. An ordinary equation of motion is second order in time, and thus requires initial conditions Zo and 
Vo for the position and velocity, but since equation (|15[) is third order in time we need an additional initial condition 
ao for the acceleration. If ao is not chosen properly we obtain unphysical runaway solutions. We can see this for the 
case of a free particle (K e xt = 0), for which the solution to (fT5"|) is 

a(i) =a Q e t/TR . (16) 

Thus, unless a — we obtain a runaway solution in which the acceleration increases exponentially in time. We can 
eliminate the runaway solutions by expressing the solution to (|15[) in terms of a Green function G(t): 

a(t) = (l/m R ) J G(t - K ext (t') dt' , (17) 

where G(t) is defined such that G(t) — > for t — » ±oo and 

G(t)-T R G(t) = 5(t). (18) 

Let G(lo) denote the Fourier transform of G(i). From equation (fl"8|) . it follows that 

G{u>)= [ G(t)e lut dt=(l + iT R Lu)-\ G{t) = {2tt)- 1 { G(uj)e- lult du J = (l/T R )e{-t)e t/TR . (19) 



We can view equation (I17p as a new equation of motion for the particle. By its construction it is a solution to equation 
(I15[) . so integrating equation (|17p subject to the initial conditions {zo,Vo} is equivalent to integrating equation (|15|) 
subject to the initial conditions {zo,vo,ao}, where ao is given by 

/>oo 

(m R r R )- 1 / e- t '/ TR K ext (t')dt'. (20) 



a 

In effect, the integro-differential equation sets the initial condition a so as to eliminate the runaway solutions. We 
can see this for the case of a free particle: when K ea;t = 0, equation ([20|) sets ao = 0. Note that since G(t) > for 
t < 0, equation (|17[) exhibits advanced effects, in which the acceleration of the particle at time t depends on the value 
of the external force at times t' > t. 

For simplicity, we have shown how the runaway solutions can be eliminated in the low-velocity limit, but the same 
method can be applied to the full relativistic Lorenz-Dirac equation (see section 6.6 of [l|, section V.6 of [3], and [a]). 
The resulting equation of motion is 

/•oo 

a»(T) = / e-V-rVm ((«(/) . a ( T ')) ^ (r ') + (m R r R )- 1 K^ xt (r')) dr' . (21) 



Note that a^{r) depends on the value of K^ xt {r') at proper times r' > t, so equation ([2T]) also exhibits advanced 
effects. 

III. EXTENDED PARTICLES IN ELECTRODYNAMICS 

We can gain some insight into the origin of the advanced effects by generalizing equation (I2ip to the case of a 
spatially extended particle of arbitrary size. For simplicity, let us consider an extended particle that is spherically 
symmetric. The current density is then given by 

J"(x) = ej(l - o(r) • (x - z(t))) /(-(as - z(r)) 2 ) w"( T ) S (v(t) ■ (x - z(t))) dr, (22) 





FIG. 1: Spacetime diagrams, which illustrate pairs of events x 11 and x'^ that give nonzero contributions to K^(t,t'). The 
thick solid lines indicate the particle trajectory, the thin solid lines indicate the planes of simultaneity for proper times r and 
t', the dashed lines indicate the light-like interval between events x^ and x ltL . (a) Conditions (1-3) are satisfied for t > t', 
corresponding to retarded effects, (b) Conditions (1-3) are satisfied for t <t', corresponding to advanced effects. 



where e is the total charge of the particle and f(r 2 ) describes the radial charge distribution (this expression for the 
current density is derived in section 7.4 of [l|). It is convenient to express the current density in the form 



J»{x) 



J»(x,T)dT, 



where 



J"(x, t) = e(l - a(r) ■ (x - z(r))) f(-(x - z(t)) 2 ) v"(t) S(v{t) ■ {x - z{t))). 



(23) 



(24) 



Note that J^(x, r) vanishes unless x^ lies in the plane of simultaneity for proper time r. We can describe the coupling 
of the particle to the electromagnetic field in terms of an action Si and Lagrangian L^. 

Si = - J A 1 " {x) J M {x) d*x = j L.dr, Li = - J A 1 " {x) J„ {x, r) d 4 x. (25) 

From the Euler-Lagrange equations, it follows that the electromagnetic force on the extended particle is 

F^(x)J v {x,T)d i x. (26) 



As before, we can separate the electromagnetic force into a component K% that describes the self-force and a component 
K? n that describes the force exerted on the particle by incoming radiation. Using equation © to solve for A${x) in 
terms of the current density, we find that the self-force is given by 

K?(t) = j F? v (x) J v (x, r) d 4 x = j Kf(r, r') dr' , (27) 

where we have defined 

K?(t,t')= [[ J a (x : T)A^ aP (x,x')jf 3 (x',T')d i xd i x', ^ a p{x,x') = {r }ct ^d> l - V %d a )D r {x-x l ). (28) 



If we assume that there is no incoming radiation, the equation of motion for the extended particle is 



(29) 



where Kp is given by equation (|27[) and K^ xt describes an arbitrary externally imposed force. This is a second-order 
integro-differential equation; it is exact, and holds for extended particles of arbitrary size. The equation of motion 
(|21[) that we derived in the previous section should approximate this equation of motion in the limit of a small particle 
(the parameters m^ = m + ms and t^ = (2/3)(e 2 /mij) that appear in equation (|2 1|) are set by the value of the 
self-energy ms, which can be determined from /(r 2 ) using equation (|A11[) ). 

In general, the self-force acting on the particle at proper time r depends on the state of the particle at proper times 
both earlier and later than r. We can understand this by examining equation (|2"8")l for K^(t,t'). From the three 
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FIG. 2: Spacetime diagrams for a uniformly accelerated particle of radius R. Black curve: particle trajectory z M . Blue curves: 
trajectories zt and z^ of the left and right edges of the particle. Red curves: trajectories L^_ and L/t. Green curves: planes of 
simultaneity for r = ±1/4. (a) aR - 1/2. (b) aR - 3/2. 



factors in the integrand, it follows that K^(t,t') vanishes unless three conditions are met: (1) there is an event x^ 
within the world-cylinder of the particle that lies on the plane of simultaneity for proper time t, (2) there is an event 
x ,fl within the world-cylinder of the particle that lies on the plane of simultaneity for proper time r', (3) events x^ 
and a/ M are light-like separated, with cc'^ earlier than x^. There are retarded effects when K^(t,t') is nonzero for 
t > t' , and advanced effects when it is nonzero for r < t' (see Figure Q}. 

We can give a concrete example that illustrates these conditions by considering a uniformly accelerated particle in 
(1 + 1) dimensions. The position, velocity, and acceleration of the particle are given by 



*"(t) 



-,- 1 D M 



e?(r), ^(r) = eS(r), a"(r)=aef(r), 



where 



Sq(t) = (cosh ar, sinh ar), 



ei(r) = (sinh ar, cosh ar). 



(30) 



(31) 



Note that eo • eo = —e.\ ■ e\ = 1 and eo • &i = 0. It is convenient to define a new coordinate system (A, u) by 

x* x (\,u)=z> J '(\)+ue?(\) = {a- 1 +u)e1(X). (32) 



The trajectory of the particle is then given by z m (t) 



'(t, 0), and the plane of simultaneity for proper time r 



is given by x^{t, u); note that the planes of simultaneity for all values of r pass through the origin at u = — 1/a. 
Equation (12~4")) for J^(x,t) generalizes naturally to (1 + 1) dimensions; from the above expressions, it follows that 



J"(x, t) = e(l + au) f(u 2 ) e^(r) 5(A - r). 



(33) 



We will assume that /(u 2 ) > for u < R and /(w 2 ) = for u > R, so the trajectories of the left and right edges of 
the particle are given by zt(X) = x M (A, —R) and z^(A) = cc M (A, +-R). 

If aR < 1, then the charge density is positive in the region between the trajectories zt and zff, and zero everywhere 
else (see Figure^). Thus, conditions (1-3) are met only if t > r'. If aR > 1, then the situation is more complicated. 
Let us define trajectories £±(A) = (A, ±|A|); these trajectories define the left and right edges of the forward and 
backward light-cones for the origin. Also, let us define i?+ to be the region between L+ and zft, and i?_ to be the 
region between z^ and Vt_ . The charge density is positive in the region R + , negative in the region i?_ , and zero 
everywhere else (see Figure HJj). Thus, conditions (1-3) are met by events in region R + if r > t', and by events in 
region R_ if r < r'. 



IV. TOY MODEL OF ELECTRODYNAMICS 



We can gain further insight into the advanced effects by comparing ordinary electrodynamics with a toy model of 
electrodynamics in which advanced effects do not occur. A complete description of the model is given in [6], but all 
the results that we will need are summarized here. 



The toy model that we will be considering describes a spatially extended particle in (1 + 1) dimensions that obeys 
Newtonian dynamics and is coupled to pair of fields E(t, x) and B(t, x), which correspond to the electric and magnetic 
fields of ordinary electrodynamics. The equations of motion for these fields are 

d t E(t,x) = 8 x B(t,x), (34) 

d t B{t,x) = d x E(t,x)-2p(t,x), (35) 

where p(t, x) is the charge density. We will assume that the charge density has the form 

p{t,x)=gf(x-z(t)), (36) 

where z(t) is the position of the particle at time t, g is its charge, and f(x) describes the charge distribution. The 
equation of motion for the particle is 

mz = F f + F ext , (37) 

where m is the particle mass, 

F/(t) = -2 / p(t, x)E(t, x) dx (38) 

is the force that the -E-field exerts on the particle, and F ext describes an arbitrary externally imposed force. Equations 
([3"4"H3"8|) give a complete description of the coupled particle-field system. Equations (j3~4"|) and (f3"5)) can be though of as 
the analogs to Maxwell's equations, and equation (j3"5)) can be thought of as the analog to the Lorentz force law. 
By analogy with electrodynamics, we can express the E and B fields in the form 

E(t,x) = E r (t,x)+E in (t,x), B(t,x)=B r (t,x)+B in (t,x), (39) 

where E r (t,x) and B r (tj_x) are the retarded fields generated by the particle, and Ei n (t,x) and Bi n (t,x) describe 
incoming radiation (see |7j). The retarded fields are given by 

E r (t,x) = d x <f> r (t,x), B r {t,x)=dt<t>r{t,X), (40) 

where 

cf> r (t, x)= D r {t -t',x- x 1 ) p(t', x') dt' dx 1 (41) 

and D r (t,x) = 9(t — \x\) is the retarded Green function for the inhomogeneous wave equation in (1 + 1) dimensions. 
Using the decompositions given in equation (|39p . we can express the force exerted by the field as Ff = Fi n + F r , where 

F m (t) = -2 / p(t,x)E in (t,x)dx, F r {t) = -2 / p(t,x)E r (t,x)dx. (42) 

Physically, Fi n describes the force exerted on the particle by incoming radiation, and F r describes the self-force 
exerted on the particle by its own retarded field. In what follows we will assume that there is no incoming radiation, 
so Fi n = and Ff = F r . 

Using equations (|40[). (|4ip . and (I42[) . we can evaluate the self- force explicitly for the case of a point particle, for 
which f(x) = 8(x): 

F r = -m 1 z{l-z 2 )- 1 , (43) 

where 7 = 2g 2 /m is a damping constant. From equations (|37[) and (|43p . we find that the equation of motion is 

z + 1 z(l-z 2 )- 1 = (l/m)F ext . (44) 



This is the toy model analog to the Lorentz-Dirac equation (J14I) for electrodynamics. If we compare equation (|44|) 
with the Lorentz-Dirac equation, we note three important differences. First, the Lorentz-Dirac equation is third order, 
but equation (|44|) is only second order, and thus does not admit runaway solutions and does not need to be replaced 
with an integro-differential equation. In the toy model the acceleration of a point particle at time t only depends 
on the value of F ext at time t, and there are no advanced effects. Second, for the Lorentz-Dirac equation there is a 
radiation damping term proportional to the time-derivative of the acceleration, while for equation (|44[) the damping 
term is a function of the particle velocity. For electrodynamics, a velocity-dependent damping term is ruled out by 
Lorentz invariance, but the toy model is neither Lorenz nor Galilean invariant; there is a preferred reference frame 
in which the equations of motion for the model are valid, and a particle moving with respect to this preferred frame 
feels a velocity-dependent drag force. Third, the mass that appears in the Lorentz-Dirac equation is the renormalized 
mass, the sum of the bare mass and the self-energy, while the mass that appears in equation (|44[) is just the bare 
mass; there is no self-energy contribution. We will explain the reason for this in section TV Al 
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FIG. 3: Spacetime diagram, which illustrates a pair of events (t,x) and (t',x') that give a nonzero contribution to F r (t,t'). 
The thick solid line indicates the particle trajectory, the thin solid lines indicate the planes of simultaneity for times t and t' , 
the dashed line indicates the light-like interval between events (t',x') and (t, x). 

V. EXTENDED PARTICLES IN THE TOY MODEL 

Let us now consider a spatially extended particle in the toy model. From equations (|4"D")) . (f4"l"T) . and (I42|) . it follows 
that the self-force for an extended particle is given by 



F r (t) = / F r (t,t')dt', (45) 

where we have defined 

F r (t,t') = -2 p(t,x)d x D r (t-t',x-x')p(t',x')dxdx'. (46) 



Equations (|45| and (|46]) are the toy model analogs of equations ((27)) and (|28f for ordinary electrodynamics. Note, 
however, that whereas K t (t,t'} can be nonzero for either t > t' or r < r', the analogous quantity F r (t,t') is only 
nonzero if t > t' . Physically, this is due to the fact that the particle obeys Newtonian dynamics, so the planes of 
simultaneity do not tilt (see Figure [31 and compare with Figure [1] for electrodynamics). In what follows, we will 
assume that the charge distribution of the extended particle is 

/(x) = (2^ 2 )- 1 / 2 e - a;2 / 2ff2 , (47) 

where a describes the particle size. If we substitute equation (|47[) into equation (|46p . we find that 

F r (t, t') = (m 7 /2V^r) 6(t ~ t') ( e -(*(«)-(0+'-0) 2 /4 CT 2 _ e -( z (t)- z (t')+t'-tf/^y (4g) 

The equation of motion for an extended particle is 

mz = F r + F ext , (49) 

where F r is given by equation (|45|) and F ext describes an externally imposed force. This is the toy model analog 
to equation (|29p . the equation of motion for an extended particle in electrodynamics. Like equation (|29p . it is a 
second-order integro-differential equation that is exact, and holds for extended particles of arbitrary size. Given the 
velocity vq at time to, together with the particle trajectory z{t) at all times t < to, equation (|49[) can be integrated to 
obtain the particle trajectory at all times. 

A. Approximate equation of motion 

We can further develop the analogy between the toy model and electrodynamics by performing a series expansion 
of the self-force for an extended particle. If we substitute equation (|4*6"]) into equation (|4"5)) and perform the integral 
over t' , we find that 

F r (t) = -2 // e(x - x') p(t, x) p(t - \x-x'\,x')dxdx'. (50) 



Let us expand p(t — \x — x f \ z x') in \x — x f \: 

oo 

p(t -\x- x'\,x') = £)(l/n!)(-l) n |z - x'\ n d?p(t,x'). (51) 

n=0 

If we assume that the particle is moving slowly, then we can neglect terms that are nonlinear in z^ n \t) = d n z(t)/dt n : 

oo 

p(t -\x- x'\,x') = p(t, x') - ^(l/n!)(-l) n \x - x'\ n z^ n \t) d x ,p{t,x'). (52) 

n=l 

Substituting equation (|52[) into equation (|50p . wc find that 

oo 

F r (t) = -m 7 ^(-l)"c„a"^")(i), (53) 

?i=0 

where v^ n ' = d n v(t)/dt n , v = z is the particle velocity, and 

en = (1/n! a") || |s - x'\ n f{x)f{x') dxdx' = ^= F( ^ ± ^ 2) (54) 



are dimensionless coefficients. The first few terms of equation (|53|) are: 

F r = —mjv — msv + mrov + ■ ■ ■ , (55) 

where ms = — (4/Y / 7r)g 2 cr and ro = —2g 2 a 2 /m. The first term is just the low- velocity limit of the self-force for a 
point particle, and the second term describes the self-energy ms- Since the self-energy is proportional to the particle 
size, it vanishes in the point particle limit, which explains why it is the bare mass, rather than the renormalized mass, 
that appears in the point particle equation of motion (1441) . 

It is instructive to compare equation (|55p with the the low-velocity limit of the electrodynamic self-force given in 
equation (TT2"|) : 

K r = —rasa. + rar^k + ■ ■ ■ , (56) 

where the dots indicate terms that vanish in the point particle limit. Note that for electrodynamics the term pro- 
portional to the velocity is not present (as we discussed before, such a term is ruled out by Lorentz invariance). 
Also, the self-energy ms and time constant To are both negative for the toy model, but positive for electrodynamics. 
These quantities also scale differently in the two theories: for the toy model, ms ~ <r and to ~ <r 2 , whereas for 
electrodynamics ms ~ 1/u and To is independent of the particle size. Thus, in the point particle limit, 

F r — > —77171;, K r — ► — m<ja + jtitqsl. (57) 

For a small enough particle, we can obtain a good approximation to equation (|53|) by truncating the series at some 
finite order N: 

N-l 

F r = -mj^2(-l) n c n cr n v ( - n) . (58) 

n=0 

Thus, we obtain an approximate equation of motion for the particle: 

AT-l 



v (l) 



+ 7 E (- 1 )" c ^ n v {n) = a/m)F ext . (59) 



n=0 



This should be a good approximation to the exact equation of motion (|49[) provided the particle is small and slowly- 
moving. Since this equation is of order N, we need initial conditions z( n )(0) for n = 0, 1, • • • , N — 1. 
We can write down the solutions to equation ([59]) for the special case F ext = 0: 



7V-1 



v(t) = J2 A k e ~ fh,t/ '> ( 6 °) 



fc=i 



where the constants A\, • ■ ■ , Ajv_i are set by the initial conditions, the constants j3i, • ■ ■ , Pn-i are the N — 1 roots 
of the polynomial 

JV-l 

p(/3)=/3-77^c„/3", (61) 

and 77 = 7cr is a dimensionless measure of the particle size. Note that some of the roots may be complex, and if so 
then the initial conditions must be chosen such that v(i) is real (the fact that there are complex roots suggests that 
there are oscillatory solutions to the equation of motion; examples of such solutions are given in Appendix IBJI . Also, 
note that if Re/3fc < and A k 7^ 0, then we obtain a runaway solution in which the velocity increases exponentially 
in time. We can write down the roots explicitly in the limit of small particle size (77 <C 1): 

fa-i=ri, (3 k = e^ k -W N - 2 Hr 1 c N -i)- l /( N -V forfc = l,--.,iv-2. (62) 

Note that for N > 4 there is at least one root that yields a runaway solution. For the case N = 3 we can calculate 
the roots exactly and write down the solutions explicitly; this is done in Appendix [Cl 

B. Integro-differential approximate equation of motion 

As for electrodynamics, we can eliminate the runaway solutions by replacing the approximate equation of motion 
([55| with an integro-differential equation. Let us define a Green function G(t) by 

JV-l 

G (D( t ) + 7 J2 (-1)™ On a n G {n \t) = 5(t). (63) 

71 = 

We will assume that N > 3, so that the polynomial p(/3) defined in equation (fuTj) can be factorized as follows: 

JV-l JV-l 

p(/3) = /? - V Yl c ^ n = -^"i II 08 - M- (64) 

Using this result, it is straightforward to show that the Fourier transform G(cu) of G(t) is 

. JV-l 

G{u) = / G(t) e iut dt = (cn^)- 1 H {iwo - /3k)- 1 . (65) 

J 1 — 1 



fe=i 



Thus, G(t) is given by 



„ JV-l 

G(t) - (27T)- 1 / G{u) e~ wt diu=Y, B >z e-> 3kt/ ' 7 0(e k t), (66) 



fc=i 
where we have defined constants B\, ■ ■ ■ , Sjv-i and ei, • • ■ , ejy-i by 



■1 ifRe/3 fc >0 



B^-lft/^-On^-A') -1 , (k ^\-l ifReS<0 (67) 

jjtk 

Using the Green function, we can express the solution to the approximate equation of motion (|59p as an integro- 
differential equation: 

r N ^ r 

v(t) = (1/m) / G(t - t') F ext (t') dt' = (1/m) J2 Bk 9( - 6k ^ ^^^ F ^ f ~ r ) dr - ( 68 ) 

•* fe=i •* 

As for electrodynamics, by expressing the solution to a higher-order equation of motion in terms of a suitable Green 
function we have eliminated the runaway solutions; for example, equation ()68j) implies that v(t) = for F ext = 0. 
Note that terms with e k = 1 only depend on the value of the external force at times t' < t, and thus describe retarded 
effects, while terms with e^ = — 1 only depend on the value of the external force at times t' > t, and thus describe 
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FIG. 4: Graph of v(t) versus t/a for an impulsive force (r) = 0.25, 
Approximate integro-differential equation of motion (|68[) with TV = 4. 



vq = 0.1). (a) Exact equation of motion (|49[) . (b) 



advanced effects. The advanced effects are not present for the exact equation of motion (|4"9f : rather, they are an 
artifact of the approximations used to obtain equation (|68[) . 

From equation (|62| . it follows that in the limit of small particle size there will be advanced effects if N > 4. Let 
us compare this with electrodynamics. Recall that the low-velocity limit of the electrodynamic self-force is given by 
equation (|56[) . It is the point particle limit of this self- force that appears in the equation of motion ()15|) , which, as we 
showed in section [Til gives rise to advanced effects. But taking the point particle limit of equation (|56|) is equivalent 
to truncating a series expansion of the self- force at order N — 3. Thus, advanced effects show up at third order for 
electrodynamics, but only at fourth order for the toy model. This can be understood by comparing equations (|55|) 
and (|S"6")) for the self-forces in the two theories. For both theories, the sign of the third-order term is determined by 
the sign of To: for electrodynamics tq is positive, corresponding to advanced effects, while for the toy model tq is 
negative, corresponding to retarded effects. Also, note that for the toy model the higher-order terms in the self-force 
vanish in the point particle limit, and thus the spurious advanced effects exhibited by equation (|68| go away in this 
limit. For electrodynamics, however, the third-order term of the self-force is independent of the particle size, and 
hence the advanced effects remain in the point particle limit. 

As an example, let us consider the special case of an impulsive force F ex t(t) — mvo6(t), and compare the evolution 
predicted by the exact equation of motion (|49[) with the evolution predicted by the approximate integro-differential 
equation of motion ((68)) . Suppose that the particle starts at rest at the origin, so z(t) = for t < 0. The moment 
after the impulsive force is applied, the velocity of the particle is v(0) — vq. We can numerically integrate the exact 
equation of motion (|4"9"|) for these initial conditions; the result is shown in Figure 2^,, where we have taken r/ = 0.25 
and vq — 0.1. As expected, there are no advanced effects: the particle does not respond to the impulsive force until 
after it has been applied. We can also describe the evolution of the particle using the approximate integro-differential 
equation (f6"5)) . If we substitute for F ext (t) , we find that 



JV-l 



"(*) = vo Yl B 



fee 



-Pkt/T 



0(e k t). 



(69) 



fc=i 



In FigureUJa this solution is shown for r\ — 0.25, vo = 0.1, N — 4. Now there are advanced effects: the particle begins 
to accelerate before the impulsive force has been applied. 



VI. CONCLUSION 



We have shown that the reason for the advanced effects in classical electrodynamics can be traced to the fact that 
the planes of simultaneity for an extended particle can tilt in such a way that the self-force exerted on the particle at 
proper time r depends on the state of the particle at proper times r' > r. We have also considered a toy model of 
electrodynamics in which the planes of simultaneity do not tilt, and have shown that it does not give rise to advanced 
effects. 

We have written down three equations of motion for electrodynamics. Equation (|29p is an exact equation of motion 
for an extended particle; it is a second-order integro-differential equation that exhibits advanced effects. Equation 
(|14p . the Lorentz-Dirac equation, approximates equation (|29[) in the limit of a small particle; it is a third order ordinary 
differential equation, and admits runaway solutions. Equation (|2ip is obtained from the Lorentz-Dirac equation by 
choosing the initial conditions so as to eliminate the runaway solutions. It also approximates equation (|29[) in the 



11 

limit of a small particle, and like equation (|29[) it is a second-order integro-differential equation that exhibits advanced 
effects. 

We have written down four equations of motion for the toy model, which can be viewed as analogs to the various 
electrodynamic equations of motion. Equation (|49[) is an exact equation of motion for an extended particle; it is a 
second-order integro-diffcrcntial equation, and does not exhibit advanced effects. Equation (|44[) is an exact equation 
of motion for a point particle; it is a second-order ordinary differential equation, and does not exhibit advanced effects. 
Equation (|59[) approximates equation (|49p in the limit of a small particle; it is an ordinary differential equation of 
order N, and admits runaway solutions for iV > 4. Equation (|68[) is obtained from equation (|59|) by choosing the 
initial conditions so as to eliminate the runaway solutions; it is a first-order integro-differential equation, and exhibits 
spurious advanced effects for N > 4. 

APPENDIX A: SELF-ENERGY OF AN EXTENDED PARTICLE 

In this Appendix we show that the self-force given in equation (|27[) can be used to obtain the correct expression for 
the self-energy of a spatially extended particle. To accomplish this, we will use equation (|27|) to calculate the self-force 
that acts on a spatially extended particle undergoing uniformly accelerated motion. For simplicity we will work to 
first order in the acceleration, and calculate the self-force at the moment at which the particle is instantaneously at 
rest. The trajectory, velocity, and acceleration of the particle are 

z»{t) = (r, ar 2 /2), «"(r) = (1, ar), o"(r) = (0, a). (Al) 

From equation (pM)) , we find that 

J"(x, t) = e(l + a • x)/(|x - z(r)| 2 ) v"(t) S(t - (1 + a • x)r), (A2) 

where t = x°. The current density is given by 

J"(x) = f J^(x,T)dT = e //(|x-z(T)| 2 )^(r)5(T-(l + a-x)- 1 t)dr, (A3) 

so 

J°(x) = p(t,x), J(x) = atp(i,x), (A4) 

where we have defined 

p(i,x) = e /(|x-a£ 2 /2| 2 ). (A5) 

Also, note that 

J»(x, 0) = (1 + a • x) p(0, x) $"„ 5{t - t). (A6) 

Substituting equation (|A6[) into equation (f27|) . we find that the self- force at r = is 



K r (0)= J p(0,x)E r (0,x)d 3 x+ f a-xp(0,x)E r (0,x)d 3 x, (A7) 

where E r (t, x), the retarded electric field, is given by 

E r (i,x) = -V^(t,x)-5 t A r (*,x), (A8) 

and A%(x), the retarded vector potential, can be obtained from the current density J^(x) via equation ((6]). After a 
lengthy but straightforward calculation, we find that 

E r (i,x) = -V f\x-x'\- l p(t,x')d 3 x'- (2/3) a f |x - x'p 1 p(t,x') d 3 x'. (A9) 

Substituting this result into equation (|A7[) , we find that the self- force is 

K r (0) = -msa, (A10) 

where the self-energy ms is given by 

m s = (e 2 /2) J J |x - xT7(|x| 2 ) /(|x'| 2 ) d 3 xd 3 x'. (All) 



It is interesting to note that the Abraham-Lorentz self- force corresponds to just the first term of equation (|A7[) . 
and gives an incorrect value of (4/3)m s for the self-energy (see @, and section 17.3 of Q). The correct relativistic 
expression for the self- force given in equation (|27p modifies this result by producing the second term of equation (|A7|) . 
which evaluates to — (l/3)ms and combines with the first term to give the correct value for the self-energy. 
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FIG. 5: Graph of v(t) versus yt for an impulsive force (vo = 0.1). Red curve: ya = 0.25. Green curve: ya = 0.5. Blue curve: 
ya = 1.0. 

APPENDIX B: EXAMPLE SOLUTIONS FOR A SPATIALLY EXTENDED PARTICLE 

In this Appendix we consider some example solutions to equation (|49[) , the exact equation of motion for an extended 
particle in the toy model. Note that one region of an extended particle can cause a change in the field that acts back 
on a different region of the particle at a later time; as we shall see, if the particle is large enough the delay between 
these events can lead to oscillatory behavior. 

As for the example in section IV B| we will assume that the particle is initially at rest at the origin, and that it 
is driven with an impulsive force F ext (t) = mvo5(t); this is equivalent to taking the initial conditions of the particle 
to be v(0) = vo, z(t) = for t < 0. We numerically integrate the equation of motion (j¥5|) subject to these initial 
conditions, and plot v(t) versus yt in FigureOU Curves are shown for three different values of the particle size a. Note 
that the renormalized mass of an extended particle is m R = m + ms = (1 — lyo j \pK)m, so for the curve with ya = 1 
the renormalized mass is negative. 



APPENDIX C: SOLUTION FOR N = 3 



Here we find the solutions to equation (|59p. the approximate equation of motion for the toy model, for the case 
./V = 3. We will assume there is no external driving force (F ext = 0), so we can express the equation of motion as 



V = —JrV + t r v, 
where tr = (m/mii)TQ, 77? = (m/mii)y, rriR = m + 7715. The solutions are given by 

v(t) = (a+ - a_) _1 ((a+v + a ) e~ a - 1 - (a^v Q + a ) e~ a+t ), 
where a± = — (1/2ta)(1 ± (1 + AynJn) 1 / 2 ), and vq and a^ are the initial velocity and acceleration. 
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